CL_H <-
function(x, s, mu, A, b, inds){
	out <- c(b - A %*% x, FB_smoother(x, s, mu, inds), mu)
	return(out)
}

